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Abstract 

We study the stationary properties of the two-dimensional pair contact pro- 
cess, a nonequilibrium lattice model exhibiting a phase transition to an ab- 
sorbing state with an infinite number of configurations. The critical probabil- 
ity and static critical exponents are determined via Monte Carlo simulations, 
as well as order-parameter moment ratios and the scaling of the initial density 
decay. The static critical properties are consistent with the directed percola- 
tion universality class. 
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I. INTRODUCTION 



Critical phenomena at absorbing-state phase transitions (i.e., between an active state and 
one in which the dynamics is frozen), are of longstanding interest in statistical physics, and 
have enjoyed renewed attention due to connections with epidemics IQ, catalytic kinetics 
surface growth , self-organized criticality 0-^] , and issues of scaling and universality P,p!0| . 
In these systems, conflict between two opposing processes (e.g., creation and annihilation), 
typically leads to a continuous transition at a critical parameter value. Such transitions 
are known to fall generically in the universality class of directed percolation (DP) ||TT|-[13|, 
although the critical behavior is modified in the presence of local parity conservation |T^-|T^. 

Another interesting case (without a conservation law), appears when the dynamics can 
become trapped in one of an infinite number (in the thermodynamic limit) of absorbing 



configurations (IN AC). Systems of this sort were introduced in catalysis modeling |[T8|-pO 
their critical properties have been studied in detail by various workers pT|-[2S[]. In one 
dimension, the pair contact process (PCP) [^], and other models with IN AC exhibit static 
critical behavior in the DP class |2^j2^, but the critical exponents associated with the 
spread of activity from a localized seed are nonuniversal, varying continuously with the 



particle density in the environment P3| , P8| , and follow a generalized hyperscaling relation 
P^ , |SD[] . The anomalous spreading can be traced to a long memory in the dynamics of the 
order parameter, p, arising from coupling to an auxiliary field that remains frozen in regions 



where p = P5| , P7| , PB | . A field theory (i.e., a stochastic partial differential equation for 
p(x, t)), incorporating this memory term reproduces the nonuniversal exponents observed 
in simulations ||3T[] . 



In two dimensions the situation is much less clear. Simulation results for a microscopic 
model with INAC ||26| conflict with studies of models exhibiting the aforementioned long 
memory [0,0 • (In particular, it seems possible that for one range of densities the spreading 
dynamics is that of dynamic percolation, while for another range there is compact growth, 
perhaps with nonuniversal exponents [0, or without well defined scaling behavior [p8|].) 
In hopes clarifying the nature of critical spreading in the presence of INAC, we propose to 
study the pair contact process in two dimensions. Our interest in the PCP is motivated by 



its simplicity, compared with the model studied in Ref. [^. The present work is devoted to 
static critical behavior, and provides the critical parameter value and the "natural density" 
in the absorbing state (defined below) needed for a detailed study of spreading. These results 
confirm that the static behavior falls in the DP class; analyses of moment ratios and the 
initial decay of the order parameter provide further support. The balance of this paper is 
devoted to defining the model and simulation algorithm (Sec. II), simulation results (Sec. 
Ill), and a brief summary (Sec. IV). 

II. MODEL 

The pair contact process (PCP) is an interacting particle system: a Markov process whose 
state space is a set of particle configurations on a lattice p2| , |33| . Each nearest-neighbor pair 



of particles has a rate p of mutual annihilation, and a rate 1 — p of attempted creation. In 
a creation attempt on the square lattice, a new particle may appear (with equal likelihood) 
at any of the six sites neighboring the pair, provided the chosen site is vacant. (Attempts to 
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place a new particle at an occupied site fail.) The PCP exhibits an active phase for p < pc, 
above this value the system falls into an absorbing configuration that typically contains a 
substantial density, 0, of particles. (Any arrangement of particles devoid of nearest-neighbor 
pairs is absorbing.) 

In our simulations, we maintain a list of the Np current nearest-neighbor pairs. At each 
step we choose a pair at random from the list, and a process (annihilation with probability p, 
creation with probability 1— p). In case of annihilation, the two particles are simply removed. 
For creation, we choose a site x at random from among the six neighbors of the pair, and 
place a new particle there if x is currently vacant. (If x is occupied the configuration remains 
the same.) The time increment associated with this step is At = 1/Np, corresponding to 
one transition per pair per unit time, in agreement with the transition rates that define the 
process. Following each change we update the list of pairs. We use a square lattice of L x L 
sites, with periodic boundaries; in the studies reported here, all sites are initially occupied. 



III. CRITICAL PROPERTIES 

To locate the critical point pc we study the size-dependence of the (quasi) stationary pair 
density p, i.e., the fraction of nearest neighbors harboring a pair of particles, in surviving 
trials, following a transient during which p{t) relaxes from its initial value of unity, p is the 
order parameter for the PCP, and as such we expect that at the critical point, 

p(p., L) ~ L-^/'^- , (1) 

while off-critical values of p should yield deviations from the power law. We studied the 
pair density in systems of size L = 10, 20, 40, 80 and 160, for times tm ranging from 10^ 
for L = 10 to 5 X 10^ for L = 160, with sample sizes ranging from 10® trials {L = 10) 
to 10^ trials {L = 160). The results (see Fig. 1) show p{p,L) following a power law for 
p = 0.2005, but clearly not for p = 0.200 or 0.201. We conclude that p^ = 0.2005(2), the 
figure in parentheses denoting the uncertainty of our estimate. The data for p = 0.2005 
yield the exponent ratio P/i'± = 0.793(5), in good agreement with the value of 0.799(2) for 



DP in 2+1 dimensions 34 



We also determined the survival probability, P{t,p,L), i.e., the probability that the 
system contains at least one nearest-neighbor pair. For finite L, this decays asymptotically 
as P ~ e"*/"^^, with the lifetime showing a power-law dependence on the system size at the 
critical point: 

rp(pe,i^)~i^"i'/"^. (2) 

The data for p = 0.2005 (see Fig. 2) yield i^\\/i'± = 1.79(1), reasonably close to the DP value, 
1.766(2). We find that p and the particle density (p also approach their stationary values 
exponentially: p{t) — Ps ^ e"*/"^ (similarly for 0), but on a much shorter time scale than that 
of P{t): T ^ Tp/lO. Analysis of the data for r at the critical point yields i^\\/i'± = 1.69(3). 

For p < Pc the process always falls into the absorbing state. The properties of this 
state are determined by the probability distribution (induced by the dynamics) on the set 
of absorbing configurations for system size L. Of interest is the particle density in the 
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absorbing state, in particular, the "natural" density, defined as the limiting value at the 
critical point: 



■'nat 



L 



lim L). (3) 



(In one dimension, it is only for this particle density that the spreading exponents take 
DP values In our simulations at pc virtually all of the trials end in an absorbing 



configuration before tm] the final particle density in the absorbing state yields an estimate 
for (j){pc,L). 

Since the dynamics of (p is tied to that of the order parameter, p, and since the excess 
particle density {4>{p) — (pnat) in a related one-dimensional model is known to be governed 
by the order-parameter exponent (3 ||2^, we expect that the leading finite-size correction to 
the particle density to be ~ L~^^'^^, just as for p. This is confirmed in Fig. 3. Linear fits 
to the data for L > 40 yield ~ 0.1480 + aL'^/"^ and p ~ bL-f^/"^, (with a = 0.9662 and 
b = 1.426), suggesting that the linear combination 0— {a/b)p will be essentially independent 
of L. This is indeed so for L > 20, as shown in the inset of Fig. 3, from which we obtain 
our final estimate, = 0.1477(1). 

Order parameter moment ratios provide another tool for assigning a model a universality 
class; in equilibrium spin systems Binder's reduced fourth cumulant has been widely used 
for this purpose A variety of ratios, involving both odd and even moments, have been 



determined for several one dimensional models with absorbing-state transitions (including 



the PGP), as well as for the two-dimensional contact process p9[. We determined the sta- 
tionary order-parameter moments mi,. ..1714 {rrij = {p-')) in order to evaluate various ratios; 
the results are listed in Table I. In Fig. 4 we plot several of the moment ratios versus L~^; 
linear fits yield the infinite-L estimates given in Table I. The latter agree quite well with the 
results for the contact process, providing further support for the PGP belonging to the DP 
universality class. Guriously, the moment ratios for the PGP appear to approach their lim- 
iting values monotonically (with L), while the two-dimensional GP exhibits a nonmonotonic 
L-dependence (see Fig. 5 of Ref. []29|). 

Finally, we analyzed the initial decay of the pair density at the critical point. In general, 
we expect the order parameter to decay as a power law, p ~ in a critical system at short 
times (i.e., t < r ~ L'^n/'^^, for which the correlation length ^ < L). Our results, plotted 
in Fig. 5, show a power-law decay with an exponent of 0.443(5). This is somewhat smaller 
than, but still consistent with, recent estimates of 6 for DP in 2+1 dimensions, which range 
from 0.4505(10) || to 0.452(1) 



IV. SUMMARY 

We studied the stationary critical properties of the pair contact process in two dimen- 
sions. On the basis of the exponent ratios and moment ratios, and the initial 
decay of the order parameter, we can assign the PGP to the directed percolation universality 
class, generic for absorbing-state transitions without a conservation law or special symmetry. 
We have noted several minor discrepancies between our results and the standard DP values, 
but expect that these are due to finite-size effects and/or a small error in p^, and do not 
reflect non-DP universality. The issue of spreading dynamics will be addressed in future 
work. 
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TABLES 



TABLE I. Ratios of order-parameter moments in the critical PCP. Entries for L = oo represent 
linear extrapolations; data for the CP from Ref. [29]. Numbers in parentheses denote uncertainties 
in the last figure. 



L 


111,2 /m\ 




771,3/ (mi 777,2) 


7714/7772 


20 


1.362(3) 


2.300(6) 


1.638(5) 


2.247(10) 


40 


1.343(3) 


2.147(7) 


1.599(5) 


2.159(9) 


80 


1.334(3) 


2.111(9) 


1.582(6) 


2.116(10) 


160 


1.327(4) 


2.086(9) 


1.571(8) 


2.093(13) 


oo 


1.323(3) 


2.067(9) 


1.56(1) 


2.07(1) 


CP 


1.326(1) 


2.080(1) 


1.569(1) 


2.093(8) 
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Figure Captions 

FIG. 1. Stationary pair density versus system size for p = 0.21, 0.202, 0.201, 0.2005, 0.200, 
and 0.195 (left to right). 

FIG. 2. Relaxation times versus system size at the critical point. Upper set: rp, the mean 
lifetime; lower set: r, associated with the relaxation of p and (p. 

FIG. 3. Particle density in absorbing configurations at the critical point versus L'^^'^-^. The 
inset is a plot of — 0.6776p. 

FIG. 4. Order-parameter moment ratios in the critical PGP. Upper set: m4/m|; middle: 
"^3/("^im2); lower: m2/m\. 

FIG. 5. Decay of the order parameter in the critical PCP. +: L = 160; •: L = 320; line: 
L = 640. 
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